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We experimentally demonstrate quantum process tomography of controlled-Z and controlled-NOT 
gates using capacitively-coupled superconducting phase qubits. These gates are realized by using 
the 1 2) state of the phase qubit. We obtain a process fidelity of 0.70 for the controlled-phase and 
0.56 for the controlled-NOT gate, with the loss of fidelity mostly due to single-qubit decoherence. 
The controlled-Z gate is also used to demonstrate a two-qubit Deutsch- Jozsa algorithm with a single 
function query. 



Quantum computation and quantum communication 
rely on excellent control of the underlying quantum sys- 
tem [H . Reasonable control has been achieved with a va- 
riety of quantum systems, with superconducting qubits 
emerging as one of the most promising candidates |2|. 
Recent experiments using superconducting architectures 
include demonstrations of quantum algorithms using two 
qubits [3] and the entanglement of three qubits y, [sj . A 
key element in these experiments is a two-qubit entan- 
gling gate, such as the V^SWAP Q and the controlled- 
Z (CZ) gates (sl, l5|. Because the CZ gate is simple to 
implement, has high fidelity, and can readily generate 
controlled-NOT (CNOT) logic [6], it likely will be an im- 
portant component in more complex algorithms such as 
quantum error correction. At present, however, the CZ 
gate functionality has only been directly tested for a sub- 
set of the possible input states. 

In this Letter, we demonstrate the operation of a CZ 
gate in a superconducting phase qubit, and fully charac- 
terize this gate as well as a CNOT gate using quantum 
process tomography (QPT). We additionally use the CZ 
gate to perform the Deutsch- Jozsa algorithm [sj, here 
with a single-shot evaluation of the function. The use of 
QPT provides a more complete gate evaluation than, for 
example, measuring the truth table for the correspond- 
ing CNOT gate [JQ, as it verifies that the gate will 
properly transform any possible input state. QPT for 
two- or three- qubit gates has been reported in NMR [9|, 
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optics |10l-ll2j. and in ion traps 
systems, QPT has been implemented for the V^SWAP 
gate with the phase qubit [15]. 

The electrical circuit for the device is shown in Fig. (TJ 
comprising two superconducting phase qubits A and B, 
coupled by a fixed capacitance C^. Each qubit is a su- 
perconducting loop interrupted by a capacitively-shunted 
Josephson junction. When biased close to the critical 
current, the junction and its parallel loop inductance pro- 
duce a non-linear potential as a function of the phase 
difference across the junction. Combined with the ki- 
netic energy originating from the shunting capacitance. 



unequally spaced quantized energy levels appear in the 
cubic potential. The two lowest levels are used for the 
qubit states |0) and |1), with a transition frequency 
{fio) that can be controlled by an external magnetic flux 
*^^x (*^ex) applied to the loop. The third energy level |2) 
is used as an auxiliary state to realize the CZ gate, as 
discussed below. 

The operation of a similar device has been reported 
previously pj. The state of each qubit is controlled 
by applying a rectangular-shaped current pulse (Z-pulse) 
or a Gaussian-shaped microwave pulse (X-, Y-pulse) to 
its bias coil. For an X- or Y-pulse, we simultaneously 
apply the derivative of the pulse to the quadrature (90° 
phase shifted) drive to reduce both unwanted excitation 
of the 1 2) state and phase error due to AC Stark ef- 
fect [13 ; the derivative scaling factor is determined from 
the nonlinearity of each qubit [l8|. This procedure en- 
ables us to use a Gaussian pulse with a full- width at half 
maximum (FWHM) of 10 ns, while maintaining accurate 
qubit control [l9| in spite of a rather weak qubit nonlin- 
earity (~ 100 MHz). Each qubit state is read out indi- 
vidually in a single-shot manner by injecting a large mag- 
nitude Z-pulse and then measuring the qubit flux with a 
superconducting quantum interference device (SQUID). 

The device was fabricated using a photolithographic 
process with Al films, AlO^^ tunnel junctions, and a-Si:H 
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FIG. 1. Circuit diagram for the experimental device, showing 
two flux-biased phase qubits coupled by a fixed capacitance 
Cc- A bias coil and readout SQUID are coupled to each qubit. 
The design parameters of the circuit are Iq — Iq — 2 /jA, 
Ca = Cb = 1 pF, La = Lb = 720 pH, and Cc = 2 fF. 
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dielectric for the shunt capacitors and wiring crossovers, 
all on a sapphire substrate. The device was mounted in 
a superconducting aluminum sample holder and cooled 
in a dilution refrigerator to ~ 25 mK. 

In the present experiment, the two qubits were biased 
so that = 7.16 GHz and /fo = 7.36 GHz when no Z- 
pulse was applied. The relaxation times (Ti) were mea- 
sured to be 510 ns and 500 ns for qubit A and B, respec- 
tively. The dephasing times determined from a Ramsey 
interference experiment (^j'^^^^^y^^ which showed Gaus- 
sian decay pro por tional to exp[—{t/T^^^^^^)'^] due to 
1/ f flux noise [20|, were 200 and 230 ns, respectively. 

Figure [2] shows the high-power spectroscopy for qubit 
B, which is used to guide formation of the CZ gate. We 
plot the escape probability of qubit B in gray scale as a 
function of the amplitude Ai of a 2 fis long Z-pulse (hor- 
izontal axis) and the frequency of a microwave X-pulse 
(vertical axis) of the same length. Both pulses were ap- 
plied simultaneously to qubit B, followed by the Z-pulse 
for the readout. In this way, we can probe the change 
in the resonance frequency as a function of detuning 
Ai. In addition to the main resonance line corresponding 
to somewhat broadened because of the large ampli- 
tude of the microwave pulse (A<l>g^ ~ 10 /i^o), a sharper 
line is observed on the low-frequency side of the main 
resonance; this corresponds to the two-photon excitation 
from the |0) to the |2) state [21]. The vertical distance 
between the main and two-photon lines is 1/2 the qubit 
nonlinearity A/ = /lo — /21, yielding A/ = 114 MHz for 
qubit A (data not shown) and 87 MHz for qubit B. 

We observe an avoided level crossing in the main 
resonance at Ai :^ 0.13 when the two qubit frequen- 
cies overlap = Here, the degeneracy of the 
\AB) = 1 10) and |01) states produces a splitting with size 
14.2±0.2 MHz, determined from a fit to the data, con- 
sistent with the designed capacitance Cc. The avoided 
crossing for the two-photon line at Ai 0.066 gives a 
splitting of 9.7±0.2 MHz, about a/2/2 times as large as 
the main resonance, as expected from a |11) and |02) in- 
teraction. The slope of the resonance between the two 
crossings is 1/2 that of /]^, as expected for the |11) state. 

Our interpretation of the spectroscopy is validated by 
a numerical calculation. Using three states for each qubit 
and the qubit design parameters, we calculate from the 



resulting 9x9 Hamiltonian [22|, l23| the energies for the 



coupled eigenstates. The energy bands, normalized to 
/o\ , are plotted versus the flux bias for qubit B in the up- 
per inset of Fig. O Here, a band is plotted only when its 
transition matrix element from the ground state is above 
a threshold, to simulate the appearance of the transition 
in the spectroscopic measurement 0. The (red) thick 
lines correspond to the |01) and |10) states, whereas the 
(blue) thin lines represent half of the excitation energy 
of the 1 11) and |02) states. The overall structure agrees 
well with the experimental data. 

As proposed theoretically by Strauch et al. [25|, the 
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FIG. 2. (Color online) High-power spectroscopy for qubit 
B. The escape probability (gray scale) is plotted versus mi- 
crowave frequency and the Z-pulse amplitude Ai. The single 
photon |0) |1) and two-photon |0) |2) transitions are 
visible, along with two avoided-level crossings. The upper 
inset shows the calculated states and eigenenergies, with the 
thick (thin) lines representing single (two) photon excitations. 
The lower inset illustrates the Z-pulse amplitude io for the 
CZ = -(SWAP) 2 operation. 



avoided crossing due to the degeneracy of the |11) and 
1 02) states can be used to construct a CZ gate, whose ac- 
tion produces no change in state except for |11) — 111)- 
By applying a non-adiabatic Z-pulse, the |11) and |02) 
states become degenerate (see lower inset of Fig. [2]). 
Initially in the |11) state, the system evolves as an 
iSWAP interaction, giving |^(t)) = cos{jAt/h)\ll) + 
ism{'yAt/h)\02), where 27 is the splitting energy of the 
avoided crossing and At the duration of the Z-pulse. Af- 
ter twice the zSWAP time At = /1/27, the system re- 
turns to the initial state |11), but with a minus sign. If 
the system starts in |00), |01), or |10), the state does 
not change since it is off-resonance with both avoided 
level crossings. A similar scheme using an adiabatic Z 
pulse has been used to successfully demonstrate a quan- 
tum algorithm [3], and the same (non-adiabatic) scheme 
has recently been used to create a three qubit entangled 
state in transmon qubits [5]. 

To experimentally determine the amplitude and length 
of the required non-adiabatic Z-pulse, we directly mea- 
sured the coherent oscillation between the |11) and |02) 
states. This (iSWAP)^ operation sequence is shown in 
Fig. [SJa): We first prepare the |11) state with a tt pulse 
to both qubits, and then apply a Z-pulse with ampli- 
tude Ai and length At to qubit B. Here only qubit B 
is probed and we adjust the measurement pulse ampli- 
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FIG. 3. (Color online) (a) Operation sequence for (b), the 
state evolution between |11) and |02) states, (b) Plot of |2) 
state probability of qubit B versus Z-pulse time At and Z- 
pulse amplitude Ai. The dashed lines correspond to the op- 
timal setting for the CZ gate, (c) Operation sequence for (d), 
demonstration of the CZ gate, (d) Plot of |1) state probabil- 
ity of qubit B as a function of icmp i'^^mp is fixed as 3 x 10~^.). 
The (red) solid and (blue) dashed curves are for qubit A ini- 
tialized to the |0) and |1) states, respectively. The vertical 
dot-dashed line indicates the value of i^mp for the CZ gate. 



tude so that the qubit is detected only when in the |2) 
(or higher) state [26]. In Fig.[3jb), we plot the tunneling 
probability P|2) as a function of At and Ai, which shows 
the expected chevron pattern. The minimum oscillation 
frequency occurs at a value of Ai that agrees with io de- 
termined in Fig.[2fa). The oscillation period to = 51.8 ns 
is also consistent with the splitting size of the avoided 
crossing. At the intersection of these two dashed lines, 
the time evolution of the state produces a minus sign, 
as required for the CZ gate. We stress that no discern- 
able increase in P|2) is observed (< 1%) at this operation 
point (At,Ai) = (to,io), confirming that we return to 
the 1 11) state after the CZ operation. 

Because the qubits themselves also accumulate phase (j) 
during the CZ pulse, the general unitary evolution from 
the gate is given by 
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By adding additional Z-pulses to both qubits, we can 
compensate these phases and even place the minus sign 
at any diagonal position in the matrix [3]. The compen- 
sation pulses are shown in Fig. [3fc), which consist of a 
fixed 10 ns pulse of variable amplitude Zcmp after the CZ 
pulse. In Fig. [3fd), we plot the tunneling probability of 
qubit B as a function of for fixed i^p. The phase of 
qubit B is measured through a Ramsey fringe experiment. 
The (red) solid and (blue) dashed curves correspond to 
qubit A being in the |0) or |1) state. They both show a 
sinusoidal dependence on i^p, but are shifted by tt from 




FIG. 4. (Color onhne) x matrices of CZ and CNOT gates, 
(a) Left panel: The real part of the experimentally obtained 
X matrix (xp) for the CZ gate, with Fp = 0.70. Right panel: 
The real part of the simulated x matrix for the CZ gate, with 
Fp = 0.67. (b) Left panel: The real part of Xp for the CNOT 
gate, with Fp = 0.56. Right panel: The real part of the 
simulated x matrix for the CNOT gate, with Fp = 0.52. The 
open boxes in the figure represent the ideal x matrix. 



each other, confirming the correct operation of the CZ 
gate. A similar experiment was done for qubit A (data 
not shown). 

The phases for the CZ gate are set by taking the values 
of icmp that give maximum probability when the control 
qubit is in the |0) state, as indicated by the vertical dash- 
dotted line in Fig.[3];d). Controlled-NOT (CNOT) gates 
are constructed by combining the CZ gate with single 
qubit rotations /7cnot = (/ Ry^'^) CZ (/ (g) Ry^^'^), 
where Ry represents the rotation of a single qubit state 
by an angle 6 about the y axis, and / is the identity 
operator. 

We evaluate the performance of these gates with QPT. 
For QPT, we prepare 16 input states in total, chosen from 
the set {|0),|l),|0) + |l),|0) + z|l)} for each qubit. After 
preparing these input states, we determine the density 
matrix of the output state with quantum state tomog- 
raphy [16], in which we measure each qubit along the 
six directions ±x, ±y and ±z of the Bloch sphere [l^. 
For each combination of QPT and QST pulses, we repeat 
the sequence 1800 times to obtain the joint qubit proba- 
bilities Pab = ^00,^0,^01 and Pn. After correcting for 
small measurement errors [l5|, we reconstruct the 16 x 16 
experimental Xe matrix from the resulting 16 density ma- 
trices With experimental noise, the Xe matrix found 
in this way is not necessarily physical, i.e. completely 
positive and trace-preserving. We thus use convex op- 
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TABLE I. Summary of performance for Deutsch-Jozsa algo- 
rithm. Deutsch-Jozsa functions are defined as fo{x) = 0, 
= 1, /2(x) = X, and fsix) — 1 — x. 



Element 




Deutsch-Jozsa function 
Constant Balanced 

/o /l /2 /s 


(OOIpIOO) + (01|p|01> 


ideal 


1 1 




measured 


0.29 0.28 0.76 0.74 


(io|p|io) + (iiHii> 


ideal 


110 




measured 


0.71 0.72 0.24 0.26 



timization to obtain the physical matrix Xp that best 
approximates Xe, as used in Ref. [15]. 

We plot the real part of Xp for the CZ and CNOT gates 
in the left panel of Figs. HJ^a) and (b). The results are 
displayed in the basis formed by the Kronecker product of 
Pauh operators {/, da^, — za^y, a^} for each qubit [6]. The 
open boxes represent the ideal x niatrix. The imaginary 
parts of Xp have very small magnitude (< 0.04 for CZ 
and < 0.03 for CNOT), and are shown in Ref. [24]. For 
both gates, we observe elements with large amplitudes 
at the proper positions. More quantitative evaluation is 
obtained by calculating the process fidelity Fp, defined by 
Fp = Tr(xiXp)7 where Xi represents an ideal x matrix. 
We obtain = 0.70 for the experimentally measured 
CZ gate, and 0.56 for the CNOT gate. For CZ gates 
with a minus sign at other positions on the diagonal, the 
measured Fp's are 0.68^.69, and 0.70 for CZqo, CZqi, 



and CZio, respectively [2j 

To understand the loss of process fidelity, we performed 
numerical simulations by solving the master equation 
using the experimental parameters [23|. The entire se- 
quence, including the QPT and QST pulses, was simu- 
lated to construct Xsim- The real part of Xsim is shown 
in Fig. m and reproduces reasonably well the reduction 
of the expected elements and the appearance of small 
unwanted elements. These imperfections are removed as 
we increase the single-qubit coherence time in the sim- 
ulation, which suggests that loss of in our system is 
mostly dominated by single qubit decoherence. We note 
that it is possible to obtain more information on the deco- 
herence mechanisms by analyzing the magnitude of par- 
ticular elements in the x niatrix [29]. 

By using these conditional gates, we can perform the 
Deutsch-Jozsa algorithm [6.] using the pulse sequence de- 
scribed in Ref. [3|. The experimental probability to ob- 
tain the correct answer is summarized in Table HI Because 
our phase qubit has single-shot readout, we can obtain 
the correct answer to a single function query more than 
70% of the time, greater than the 50% probability for 
a classical query and guess. We stress that no calibra- 
tion for the measurement error is applied here. The full 
density matrix of the final state is given in Ref. [24 1. 

In conclusion, we have demonstrated CZ and CNOT 



gates in capacitively-coupled phase qubits using the 
higher-energy |2) state. Quantum process tomography 
measures a x niatrix that is in good accord with predic- 
tions, which is a definitive test of proper gate operation 
for any input state. 
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I. CALCULATION OF THE ENERGY BANDS AND TRANSITION MATRIX ELEMENTS 

We calculated the energy band of the capacitively-coupled flux-biased phase qubits by diagonalizing the following 
9x9 Hamiltonian [l|,[2|, 



CN ■ /o \ /o 




- . where f-^-' (/f^^) is the flux-dependent transition frequency between ith and jth state of the qubit A (B), and g is the 
^ ' coupling energy between the qubits. The last term in the Hamiltonian is based on a^a^-type coupling of the two qubits. 
To calculate the transition matrix from the ground state \g) to the excited state |e) in the spectroscopy experiment, we 

{e\a^ -^a\i){i\a^ a\g) 2 



calculated the transition matrix element of \{e\a^ -\- a\g)\'^ for one photon excitation and ^ 



d • ' ~^ Ee — Ei — hfd 

for two-photon excitation [3], where a (a^) is an annihilation (creation) operator for the harmonic oscillator, Ei is 
the energy gap of the state \i) from the ground state, and /d is the frequency of the /i-wave drive, which is set to be 
Ee/2h in the calculation. 



II. X MATRIX FOR ALL GATES 



^ ■ In Fig. [HH the physical x matrix Xp is plotted for all the CZ and CNOT gates. 
00 

O ■ III. THE DIFFERENCE BETWEEN Xp AND Xe 

in 

' We checked the difference between Xe (the experimental x niatirix) and Xp (the physical x niatirix) by histogram- 

. ming the differences in the peak height A = Xe — Xp of each of the 256 matrix elements in the real part [4]. We fit 

*^ ■ it by Gaussian aexp(-AVcr^) as shown in Fig. [H2 The obtained cr are 0.0020 for CPn and 0.0017 for CNOT gate. 



which implies Xe and Xp are close. 



IV. SIMULATION OF QPT 



To simulate QPT, we solved the standard master equation, p = —{i/?i)\H^p\^C[p\^ where ?^ is a 9 by 9 Hamiltonian 
for capacitively coupled phase qubits under rotating wave approximation, and C[p] = ^ ^ ^)P^] ^ ~ ^ ^)P~ 

i=A,B j = l,2 

-p>C^- ^>C* . Here, for example, Ci = ua/ \/T^ and C2 = <^a^a \/2/T2 describe the relaxation and dephasing for qubit 

A, respectively [5|. Experimental Ramsey interference shows Gaussian decay, which is not reproduced by the above 
master equation. Thus, in order to approximate this situation, we used an effective T2 that depends on the length 

of the control sequence for a particular experiment tgeq- In particular, we used T2 = T^^"^^^^ jt^^^ in the simulation 
in order for both the Gaussian decay and exponential decay to give the same decay factor at tgeq. The actual tgeq is 
101.8 ns in QPT for CZ gates and 141.8 ns for CNOT gate. In Fig.[S3l we show the simulated x matrix of all CZ and 
CNOT gates. 
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FIG. SI: Xp of (a) CZoo (b) CZoi (c) CZio (d) CZ = CZn and (e) CNOT. Process fidelity Fp are 0.68, 0.69, 0.70, 0.70 and 
0.56, respectively. 
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FIG. S2: Histogram of the differences in the peak height of each of the 256 matrix elements in the real part of x matrix. Data 
is for (a) CZ and (b) CNOT. The solid curves are a Gaussian fit to the data. 



V. EXPERIMENT ON DEUTSCH-JOZSA ALGORITHM 



Figure [S4lf a) shows the pulse sequence for the Deutsch-Jozsa algorithm. The four different two-qubit gates Ui 
correspond to the four Deutsch-Jozsa functions, which we want to determine by a single quantum evaluation of the 
function. They are given by 

Ui=I^Rl , 

t/2 = (/ rV^ri) cz^^ {1 i?^/') , 

The sequence is same as that used in Ref. [H except that Ry^'^ pulse was not applied to qubit B before the tomography. 
This makes the final state a superposition state. Also, in order to shorten the total sequence time, the last R^^'^ on 
qubit A was applied before the Ui part finishes. The real part of the final density matrices are plotted in Figs [S4F b)-(e). 
No calibration for the measurement error is applied here. 
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FIG. S3: Simulated x matrix of (a) CZoo (b) CZoi (c) CZio (d) CZ and (e) CNOT. Process fidelity Fp are 0.67, 0.67, 0.66, 
0.67 and 0.52, respectively. 
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FIG. S4: (a) Pulse sequence for the Deutsch-Jozsa algorithm, (b)-(d) Real part of the density matrix of the final state for four 
Deutsch-Jozsa functions. Open boxes represent the ideal density matrix. 



